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Abstract 

It is a well known analytic result in general relativity that the 2-dimensional area of the apparent 
horizon of a black hole remains invariant regardless of the motion of the observer, and in fact is 
independent of the t = constant slice, which can be quite arbitrary in general relativity. Nonethe- 
less the explicit computation of horizon area is often substantially more difficult in some frames 
(complicated by the coordinate form of the metric), than in other frames. Here we give an explicit 
demonstration for very restricted metric forms of (Schwarzschild and Kerr) vacuum black holes. In 
the Kerr-Schild coordinate expression for these spacetimes they have an explicit Lorentz-invariant 
form. We consider boosted versions with the black hole moving through the coordinate system. 
Since these are stationary black hole spacetimes, the apparent horizons are two dimensional cross 
sections of their event horizons, so we compute the areas of apparent horizons in the boosted space 
with (boosted) t = constant, and obtain the same result as in the unboosted case. Note that 
while the invariance of area is generic, we deal only with black holes in the Kerr-Schild form, and 
consider only one particularly simple change of slicing which amounts to a boost. Even with these 
restrictions we find that the results illuminate the physics of the horizon as a null surface and 
provide a useful pedagogical tool. As far as we can determine, this is the first explicit calculation 
of this type demonstrating the area invariance of horizons. Further, these calculations are directly 
relevant to transformations that arise in computational representation of moving black holes. We 
present an application of this result to initial data for boosted black holes. 
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I. INTRODUCTION 



Apparent horizons (AH) were first introduced by Penrose and Hawking [1], [2]. An 
apparent horizon is defined as the outermost marginally trapped surface on a given (partial) 
Cauchy slice. It is a topologically spherical 2-dimensional surface on which the expansion of 
the outgoing null rays orthogonal to the surface is zero [3] . Thus, it is a surface where gravity 
is so strong that putative outgoing null rays can only "hover" against the gravitational force. 
Unlike event horizons, which are globally defined as the boundary in spacetime between 
null geodesies that escape to infinity, and those that fall into the singularity, apparent 
horizons are local objects, computable at one instant of time, hence much more accessible 
in numerical simulations. The locations of event and apparent horizons coincide only in 
stationary spacetimes. In the stationary Kerr spacetime in Boyer-Lindquist coordinates, the 
horizon is located at radial coordinate r = r+ = M+y/M 2 - a 2 , where M is the mass of the 
black hole and a is the spin parameter for the Kerr black hole given by a = J/M, with J 
being the angular momentum of the black hole; for a = we have the static Schwarzschild 
black hole. Here and henceforth, we use Newton's constant G = 1, and speed of light c = 1. 

In this paper we consider only Kerr spacetimes in Kerr-Schild (KS) coordinates as given 
in Eq. (1) below. This form of the metric contains a "natural" Minkowski background, 
and hence a natural definition of a Lorentz boost [4]. It is found (cf. [5], [6], [7]) that the 
apparent horizon of a black hole will appear distorted in these coordinates when boosted; the 
longitudinal coordinate direction undergoes a Lorentz contraction. However, this is an effect 
only in coordinates; the point of this paper is an explicit calculation to show that the area 
of the apparent horizon 2-surface, recomputed in the spatial frame of the boosted observer, 
remains unchanged, that is: Area = 4tt (r+ + a 2 ) for the Kerr case and Area = 16tiM 2 for 
the Schwarzschild black hole. This result is of course necessary on general principles. 

The invariance of the area depends on the observation that the event horizon of a station- 
ary black hole is a null 3-dimensional submanifold of the spacetime with vanishing expansion. 
And null surfaces naturally remain null under Lorentz transformations. In fact, the area of 
any 2-dimensional cross section of the horizon remains invariant under any redefinition of 
the 3-space t = constant (that is legitimately spacelike). Two cross sections of the event 
horizon that differ by a redefinition of t = constant slice can be put in a pointwise 1-to-l 
correspondence along the null generators of the horizon. These null offsets do not contribute 
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to the area which is transverse to the null generators. We give a quick derivation of the 
Schwarzschild situation and then present the most general calculation for these spacetimes, 
namely, the Kerr black hole boosted along an arbitrary direction. 

The Kerr vacuum solution to Einstein's equation can be written in a special form called 
the Kerr-Schild form of the metric. This form is, in general ([8], [9], [10], [11]), 

= f]p,v + iHl^lv (l) 

where if is a function of spacetime coordinates, T}^ is the Minkowski metric of flat spacetime 
and P is a null vector with respect to both g^ v and r}^. Clearly, this is a special form, and 
the metric of a general spacetime cannot be put in this form. But the Kerr vacuum black 
hole can be so written. Under a Lorentz boost (a coordinate transformation with the form 
of a Lorentz transformation on the t, x, y, z coordinates describing the flat space with metric 
rj^dx^dx"), the Kerr-Schild metric will preserve the general form that it has in Eq. (1). 
We will place overbars on coordinates in the unboosted frame. In section III, we will show 
the area invariance for a boosted Kerr black hole by performing a coordinate transformation 
to facilitate boosting the spacetime, followed by another coordinate transformation that 
simplifies extracting the 2-dimensional metric by restricting to the horizon. With the 2- 
dimensional metric we straightforwardly compute the horizon area. 

The special case of the nonspinning Schwarzschild (i.e. spherical) black hole provides an 
illuminating guide to the features of the full Kerr case. Eq. (1) for this case is 

gaudx^dx v = —di 2 + dx 2 + dy 2 + dz 2 H — {di + dx + dy + dz) 2 (2) 

r 

which, in cylindrical coordinates (f\\,fx, <f) C yi) can be written as 

g^dx^dx" = -dt 2 + df 2 + dx\ + dy 2 ± + — ( di + ydf l{ + ^ydx± + ^dy± J . (3) 

where x±_ = f±cos (j) cy i and y± = f±sin (f> cy i. The coordinate system (fn, f±, 4>cyi) aligns 
f|| with the axis of the cylinder parallel to the boost direction (3; fj_,4> C yi are the polar 
coordinates of the circular plane orthogonal to the axis of the cylinder (see Fig. (3)). Note 
that x 2 + y 2 + z 2 = fjj + fj_. 

The boosted (unbarred) coordinates are related to the unboosted frame by 

t = 7 {t-l3r\\) 
f|| = 7(rj| —fit) 

r± = r ± , <j) cy i = 4> cyl . (4) 



The boost parameter is (3 = v/c, and 7 = (1 — /3 2 )~ 1 / 2 ; both are defined as usual in 
the background Minkowski spacetime. The apparent horizon is defined in a given 3-space 
(t = constant) and the horizon area will be independent of t, so we take t = 0. The t = 
(boosted) 3-metric is 

ds 2 \ t=0 = dr 2 + dx 2 ± + dy\ + ~^-{— ry/3dr\\ + f\\df\\ + x±dx± + y_i_dy±) 2 . (5) 

We have strategically kept some terms expressed using unboosted (barred) forms. They can 
be straightforwardly substituted using Eq. (4). In this form, however, we can easily restrict 
the metric to the horizon surface, since 

fdf = f\\df\\ + x±dx± + y±dyj_ = (6) 

on the horizon where f is a constant (= 2M). Thus on the horizon: 

ds 2 



\t=0, f=2M 



2M 

f 

(dfl + dx\ + dy]_) \ f=2M ■ (7) 



7" 2 dfjj + dx\ + df ± + ^-{(3 2 df 2 ) 

1^2 , J-2 , J„-;2 



In Cartesian (x, y, z) coordinates this would look like 

ds 2 \ t=0 f=2M = (dx 2 + dy 2 + dz 2 ) | f=2 M ■ (8) 
This can be put in a more familiar form using spherical coordinates (f, 9, <p) which now gives 

ds\ =Q) - =2M = (2M) 2 (d9 2 + sin 2 6d^ 2 ) . (9) 

Thus the area of the horizon is 4ir(2M) 2 as expected. Importantly, note that Eq. (7) 
describes the boosted apparent horizon; the simple form (Eq. (8)) that allows immediate 
evaluation of the surface area is the expression of this area in terms of coordinates appro- 
priate first of all to the unboosted frame. On the horizon the contribution from the time 
transformation exactly cancels the Lorentz contraction of fu. 



II. NUMERICAL RESULTS 



Before looking at the horizon of a boosted spinning black hole, we demonstrate some 
numerical applications of these concepts, concentrating in this section on only nonspinning 
black holes. Recent breakthroughs in numerical relativity ([12], [13], [14], [15], [16], [17]) have 
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enabled the community to investigate various physical scenarios involving interacting black 
holes. There are many different approaches to numerically evolving the physical system. The 
use of a particular structure, puncture initial data ([18]) has become ubiquitous for numerical 
codes. Puncture initial data are conformally flat. Solution of the constraint equations 
(elliptic equations describing a nonlinear generalization of Newtonian gravity) produces a 
mathematically correct configuration. But if boosted, the puncture is not physically relaxed, 
so when the solved (mathematcally correct) data are evolved, the black hole emits short 
wavelength gravitational radiation. Some of this spurious radiation propagates out to infinity 
and some falls onto the black hole, increasing the horizon mass. 

One can instead use superposed Kerr-Schild ([19]) initial data. This takes the Kerr-Schild 
metric for a single black hole and creates a background metric for two black holes by adding 
a second 'mass term' to the flat background: 

g.u^v^ + m^ + m^. (io) 

Here Hi, H 2 are scalar functions that depend on coordinates from the centers of each black 
hole as well as the black holes' masses and spins. They are identical in form to single 
black hole terms centered at the locations of the two holes (cf. Eqs. (1), (2), and also Eq. 
(11) below; there is also a prescription for superposing the momentum associated with this 
combination, in the initial data). Although Kerr-Schild initial data exactly solve Einstein's 
equation for a single boosted black hole and thus satisfy the constraint equations, this is 
not the case for superposed Kerr-Schild, which is only an educated guess. However, by 
starting out with this initial guess as a conformal background metric (in the same sense that 
puncture data has a flat conformal background), one can solve the constraint equations, so 
Kerr-Schild data can be adjusted to become proper initial data. The solution of the elliptic 
initial data equations modifies the configuration to be an exact (modulo numerical error) 
description of a gravitational configuration. In practice, unless the black holes are very close 
together, the correction for superposed Kerr-Schild data is small; less than one percent. 

The code being developed at University of Texas Austin is called openGR [20]. Among 
the suite of programs comprising openGR, there is a finite element initial data code, which 
can produce either puncture or superposed Kerr-Schild initial data. The evolution code 
treats the dynamics of binary black hole systems and the extraction of gravitational waves 
from the merger of the black holes. The code is a fourth order accurate adaptive mesh 
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refinement code with sixth order interpolation between coordinate patches. 

The total mass/energy of the spacetime is given by the ADM mass Madm ([28]) computed 
at spatial infinity (numerically, "near" the grid boundary). The ADM mass corresponds to 
the apparent Newtonian mass measured at large distances from the sources, measured for 
instance by observing the period of distant satellites around the central mass. Suppose the 
individual black hole masses are given by Kerr-Schild mass parameters mi and m-i. Then 
the background gives an ADM mass Madm bkgd = Trii +rri2. As noted, solving the constraint 
equation changes the superposed Kerr-Schild data slightly, so the solved ADM mass closely 
approximates Madm + rn 2 , though it does have some dependence on the parameters 

of the data, particularly on the separation of the black holes. 

Of interest in the design of data is the binding energy of the configuration. We can 
compute this as the measured ADM mass minus the intrinsic mass of the constituent black 
holes. The difficulty lies in defining an intrinsic black hole mass. We choose the horizon 
mass. (For nonspinning black holes, we have Mh = (Ah/IQtt) 1 ^ 2 , where Ah is the area of 
the apparent horizon; openGR includes an apparent horizon finder.) Classically the area of 
the horizon can increase, but we also know that the horizon area is an adiabatic invariant; 
it is only slightly affected by slow motions. "Slow" means slow compared to the normal 
frequencies of oscillation of the hole, which are high frequency; the lowest frequency is on the 
order / ~ (20M#) -1 , and most frequencies in binary evolution are lower than this frequency. 
Hence we are confident that the apparent horizon provides an (almost) constant intrinsic 
mass. Binding is indicated by = Madm — (Mm + M^i) < 0. If the data describe black 
holes in motion, then the kinetic energy also contributes (positively) to the total energy. 
For a boosted black hole the ADM mass acquires a factor 7: Madm — > jMadm where 
Madm is m, the metric mass parameter in the single hole case. Thus we expect that 
for a given boost parameter 7, the binding energy may be negative (i.e. bound) if the 
holes are close together, but positive (unbound) if the data are set with the black holes far 
apart. Furthermore, for the nonlinear small separation limit (and/or for significant 7) cases, 
Newtonian arguments become obscure because of the change in metric due to the presence 
of the second hole, and due to coordinate ambiguities. 

We construct an equal mass binary black hole system (nonspinning Schwarzschild black 
holes) with initial coordinate separation r. The configuration is axisymmetric; the black 
holes are boosted toward or away from each other with Lorentz boost velocity (3 = v/c (or 
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instantaneously at rest with (3 = 0); the boosts are equal but opposite in the computational 
frame. The axisymmetry allows extremely high resolution computational simulation. The 
code is a finite element code, with an adaptive resolution of 1/100 Madm near the holes and 
1 Madm at the outer boundary. The computational domain is a sphere of radius 256 Madm- 

We plot the the negative of binding energy — Ef, of the binary in figure 1. Here, we define 
the binding energy to be E\, = Madm — 2Mh- Since the configuration is axisymmetric, the 
black holes have the same horizon mass Mh, hence the factor of 2. We display our results 
in units of the total parameter ADM mass which is normalized to equal 1 (Madm bkgd — 
mi + ni2 = 0.5 + 0.5 = 1). Here r is the coordinate separation between the two black holes 
also given in units of Madm (e.g. r — 10 translates to r = 10 G Madm/c 2 ). The binding 
energy scales as 1/r at the Newtonian limit; this Newtonian limit is plotted as a red straight 
line in Fig. 1. Bonning et al. [21] had analytically predicted this Newtonian limit (see that 
paper for details). Previous computational work by Hawley et al. [22] failed to show the 
Newtonian limit, because of insufficient domain size to eliminate outer boundary effects. We 
clearly see that for every rest configuration the binding energy for large separations agrees 
with the Newtonian prediction ([21]), but there is a deviation to stronger binding for closer 
coordinate separation. We will have more to say about this in a future paper. The cause for 
this will be discussed below as it related to the distortion of black holes' horizons near each 
other. One can in principle use expressions from post-Newtonian theory to give the next 
order correction to Ef,. These terms scale as ((Mass)/r) 2 . We have begun studying these 
higher order corrections. 

It is of interest to understand how the binding energy is achieved in the initial data. 
Fig. 2 is a plot of Madm and horizon mass M H versus 1/r for boosts of (3 — 0, 0.1, 0.5 
represented by the red, green and black curves for Mh , respectively. For the Madm versus 
1/r plot, we use a blue solid line, red "x" marks and pink dashed line for f3 = 0, 0.1, 0.5, 
respectively. Note the confirmation of the analytical expectation above that the ADM mass 
is essentially constant for the binary pair regardless of the coordinate distance between them. 
However, although we construct all data with the same parameter values m, we see different 
constant ADM masses for different \/3\ (motion with the same \j3\ together or apart yields 
the same ADM mass, constant across the possible separations). This is because the ADM 
mass scales as 7 Madmo for a boosted black hole. Thus, for example, the ratio of ADM 
masses between the pink dashed line (j3 = 0.5) and the blue line (f3 = 0) in Fig. 2 should 
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FIG. 1: Negative of the binding Energy — Ej, versus the inverse coordinate separation 1/r for 
the cases with boosts speed (3 = 0, 0.1, 0.5 represented by the green, black and blue curves, 
respectively. The red line is the Newtonian binding energy which scales as 1/r. Ideally, it should 
be tangent to the (3 = curve (green) at large r (1/r — > 0) but here it is slightly shifted due 
to numerical errors. As can be seen in the figure, the binding energy matches the Newtonian 
limit very well for large separations (1/r — > 0), it grows faster than 1/r as the black holes are 
closer (1/r — > oo). This is due to changes in horizon masses because of the distortions induced 
by the black holes on each other. It (—Eb) also becomes more negative for large boosts reflecting 
the unbound nature of distant rapidly moving black holes. The kinetic energy of the black holes 
overwhelms the negative potential energy. 

be (1 - 0.5 2 )~ 1/2 = 1.154. This is easily seen in Fig. 2. We estimate the numerical error of 
about one percent in this quantity by looking at the ADM mass for the (3 = case (blue 
line) which, in principle, should give Madm = 1 but actually is located slightly higher at 
M ADM = 1.01. 

Though Madm stays almost constant for differing separation, the binding becomes 
stronger for smaller separation, even in the Newtonian limit, of course. As described in 
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FIG. 2: Horizon mass Mh and ADM mass Madm versus inverse distance 1/r for boost speeds of 
j3 = 0,0.1,0.5. Madm (the upper, approximately parallel curves) is given by the the pink dashed 
line (v=0.5) and by the closely overlapping blue line (v=0) and the red tick marks ("x"). As 
expected, the ADM mass remains constant regardless of the separation r but varies as 7 Madm 
for varying boost speeds (3. Mh (the lower curves) is represented by closely overlapping red and 
green curves for the v=0 and v= 0.1 cases, and by the higher black curve for v= 0.5. Note that 
the horizon mass grows larger as the black holes are nearer i.e. as 1/r ^ 00. The horizon mass 
is invariant under boosts. For r > 10 (i.e. 1/r < 0.1) the horizon mass curves for different boosts 
overlap perfectly. Apparently because of the nonlinear interaction of the black hole geometries in 
the full solution, for larger boosts and for small separations the horizon mass does increase slightly. 

[21], when the parameter m is held constant for each hole, the horizon area of the con- 
stituent black holes increases with decreasing separation. The modification of the geometry 
by the other black hole modifies the horizon area so that it is no longer the 16irm 2 which 
would be computed for an isolated hole, but 16nM H with Mh 7^ m. If we imagine the 
initial data constructed by adiabatically moving the holes from infinite separation, it would 
be this mass Mh which is adiabatically invariant. This was predicted analytically by [21] 
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for nonspinning, instantaneously nonmoving black holes; it was predicted qualitatively for 
moving black holes. 

The horizon mass is expected to remain invariant under boosts, and in single boosted 
black holes this is what we observe. But for fully solved data - the result of solving a 
nonlinear elliptic system, Fig. 2 shows that the horizon mass in the f3 = 0.5 case is somewhat 
above that of the (3 = one for close separations ('close' meaning black hole separations less 
than r = 10). Indeed, for r > 10 (i.e. 1/r < 0.1), the overlap of the red, green and black 
curves is perfect to within less than one percent error. This is an interesting result depending 
on both the boost and the separation. The growth of the horizon area for large boosts is 
an effect due to the proximity of the two black holes. For sufficiently large separations, the 
boost does not change the horizon mass, hence the horizon area. 

III. BOOSTED KERR BLACK HOLE 

We return to the analytic study of black hole horizons, now including spin. Strong as- 
trophysical evidence supports the existence of spinning (Kerr) black holes ([23], [24], [25]); 
manipulating description of this spacetime is a frequent task in computational astrophysics. 
The angular momentum of the spinning black hole automatically selects a preferred direc- 
tion and the Kerr hole is axially symmetric around the spin axis. Written in Kerr-Schild 
coordinates the Kerr spacetime formally admits a boost. 

We will begin with the unboosted Kerr metric written in standard Kerr-Schild coordi- 
nates. We will then rewrite the metric in cylindrical coordinates where the symmetry axis 
of the cylinder points toward the boost direction. (We transform to cylindrical coordinates 
only to facilitate the boosting of the spacetime.) Once the spacetime is boosted, we will 
look at the spatial 3-metric on a (boosted) t = constant hypersurface. Since we are ulti- 
mately interested in the 2-metric we will perform one final coordinate transformation from 
cylindrical to spheroidal coordinates and consider f = r + (the expected horizon location). 
Once we have our 2-metric, we will compute the area of the apparent horizon and show that 
it indeed equals the unboosted, stationary value, which is Area = Air (r 2 ^ + a 2 ) . 

The Kerr spacetime in Kerr-Schild coordinates is ([8], [9], [10], [11]): 



ds 



2 



-dt 2 + dx 2 + dy 2 + dz 2 + 



2Mf 3 



dt + 



rx + ay 




r 2 + a 2 
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-dt 2 + dx' z + dy A + dt 



2Mf 3 



r , , , . a(ydx — xdy) zdz 
dt + — - — {xdx + ydy) + yy yj + -34 

r Z _|_ r 2 _|_ r 



r 4 + a 2 z 2 

where we rewrote the Kerr metric in the second line in a form that will be useful in the 
following. In the a — > limit, we recover the Schwarzschild metric in Kerr-Schild coordinates. 
The radial coordinate f is related to the fundamental coordinates x, y, z by the equation of 
an oblate ellipsoid 

* 2 + y 2 + t = 1 ( i2) 

r 2 + a 2 r 2 
which is equivalent to a quadratic equation in f 2 

f 4 -f 2 (x 2 + f + z 2 -a 2 ) -a 2 z 2 = , (13) 



and the horizon is located at r = r + = M + \/M 2 — a 2 . Eq. (12) motivates spheroidal 
coordinates: 



.r 



\fr 2 + a? sin cos i 



y = Vr 2 + a 2 sin 6* sin (14) 
z = fcos#. 

We now explicitly reintroduce the cylindrical coordinates (f\\,f±, 4> C yi) of the previous section 



x — f || sin Op cos <f>p + f± (cos Op cos <f>p cos (j) cy i — sin <j)p sin 0^) 
y = f 1 1 sin 6p sin 0/3 + fj_ (cos 6p sin 0/3 cos 4> cy i + cos 0/3 sin C ^) 
z = cos Op — f ±sin Op cos 4> cy i 



(15) 



The angles Op, (ftp specify the direction of the Lorentz boost f3 in spherical coordinates based 
on x,y,z: (3 = (/3sin Op cos 0/3, /3sin ^ sin (f>p, (5 cos Op). With the coordinate transformation 
in Eq. (15) the Kerr metric becomes 



ds 2 



-dt 2 + dr 2 + dr\ + r 2 ± > 
2Mr 



% y i 



r 1 



+ a 2 (r 1 1 cos Op - r± sin Op cos <p cy i ) ' 



x 



dt 



r 2 + a 



_ zdz a (smOp [r±dru - ru d(r± sin0 CJ// )l - cos 0pr\d(f) cy i) 

[xdx + ydy) + — + ^- — 

a r r z + or 



(16) 
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FIG. 3: The tilted cylindrical coordinates (vu,f±,</>cyi) along with the radial coordinate r and 
Kerr-Schild Cartesian coordinates (x,y,z). The vector (3 points along the boost direction, which 
is parallel to the symmetry axis of the cylinder. 

We now carry out a boost along the selected cylindrical axis. Unbarred coordinates will 
denote the boosted observer frame. They are related to the barred rest-frame coordinates 
via Eq. (4). After boosting this metric, we will look at it on an arbitrary t = constant 
hypersurface, which we take as t = since this choice simplifies the expressions, to project 
out the spatial geometry of the hypersurface in which the apparent horizon lies. This will 
leave us with dt = —^f3dr\\ and f\\ = 77*11 . With these changes substituted into Eq. (16) we 
obtain the spatial part of the boosted Kerr metric on a t = hypersurface: 

ds 2 \ t=Q = dr\ + dr\ + r 2 ± d<j) 2 cyl 

2Mf 

+ 2 

r 4 + a 2 (77*11 cos 9/3 - r± sin dp cos <j) cy i) 
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X 



r 2 + a 



,_ _ ,_s zdz a (sin ^7 [ridrn - r\\ d(r± sin</w 
2 (xax + ydy) H — = — I 



r 2 + a 2 



where a few terms involving x , y , z , r in Eqs. (16), (17) were left untouched with the next 
step in mind. If one wishes, one could also write all of these terms as functions of ru, r± and 



Let us remember what we are after; the 2-metric of the boosted geometry projected out 
by the condition f = r + . Eq. (12) implies 



xdx + ydy zdz 



x 2 + y 2 
'f 2 + a 2 ) 2 



z 2 



f 1 



— rdr 



at 



(18) 



Thus if df = (e.g. on the horizon, f = r+), the left hand side of Eq. (18) vanishes. This is 
the analogue of Eq. (6) for the Schwarzschild case of Section I. This simplification reduces 
the complexities of Eq. (17) substantially: 



ds 2 



t=0, r=r+ 



dr?, + dr 2 , 



r 2 , d(t>l,, + 



2Mr. 



J cyl 



r++a (7^|| 



-r x sin 0,3 cos <p cyl 



X 



—j(3dr\\ + 



afsinOp"/ [rxdr||— r|| d(r ± sin <f> cy i)] -cos 6pr^_d<f> cv i 



(19) 



r=rj r 



However, it is difficult to translate the horizon condition f = r + into something meaningful 
in cylindrical coordinates. Therefore, we must rewrite Eq. (19) in spheroidal coordinates to 
impose the condition f = r + to extract the 2-metric of the apparent horizon. We do this 
by going back to Eqs. (15) and rewriting them as a matrix equation for both boosted and 
unboosted coordinates 



y 



( r± cos <p cy i ' 



M 



r ± sin (p cy i 



( r ± cos 4> cy i ' 



M 



r ± sin (p cy i 
yr\\ 



(20) 



where the components of the matrix M can be determined from Eqs. (15). The radial 
coordinate f in Eq. (14) is related to the Cartesian and cylindrical coordinates via x 2 + 
j_ + f jj = f 2 + a 2 sin 2 6. Since f\\ = jr\\ on the t — hypersurface, we also have 



y 2 + z 2 



r\ + 7 2 r 2 = r 2 + a 2 sin 9 (cf. [8]). Setting Eq. (20) equal to Eq. (14) and multiplying by 
M -1 , we get 

^ y/f 2 + a 2 sin 6 cos <f> 



( r± cos 4> cy i ^ 



r ± sin C2/Z 

7T|| 



M" 1 



V / r 2_ +~a 2 sin 6* sin • 
fcos8 
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Expanding this we obtain 



r± cos <f>cyi = cos Op cos <$>p\J f 2 + a 2 sin 9 cos + cos 9p sin (fipVr 2 + o 2 sin # sin — sin 9@ f cos 9 
rj_sin0 c ^ = — sin (ppV f 2 + a 2 sin 9 cos + cos 0/jV f 2 + a 2 sin 9 sin (21) 
7r|| = sin 9p cos (ppy/f 2 + a 2 sin # cos + sin 9p siiKppy/f 2 + a 2 sin 9 sin + cos#g fcos9 . 

In the limit 9@ — 4>p = 0, the equations above reduce to Eq. (14) with the cylindrical 
coordinates replacing (x, y, z). This is the case of boosting along the z-axis, and we briefly 
treat that here before proceeding. For boost along the z-axis, with f — r+ (i.e. on the 
horizon) we have 



lt=0, r=r+ 



dx 2 + dy 2 + dz 2 , 



ri (ri + a 2 ) 



ri + 7 2 a 2 z 2 



-'jpdz + 



a(ydx — xdy) 
ri + a 2 



1 2 



(22) 



r=r+ 



Because of the boost in the z-direction, only the terms involving z (—yfidz in the numerator 
and 7 2 a 2 z 2 in the denominator) differ from the unboosted case. In Eq. (22) we still have to 
evaluate some of the terms at f = r+. Using spheroidal coordinates 



\l r\ + a 2 sin 9 cos < 
y = \l 'ri + a 2 sin sin ( 



(23) 



•yz = r + cos 9 



we get 



t=0, r=r_|_ 



(r 2 +a 2 )(cos 2 #d# 2 + sin 2 



+ -± sin 2 #d# 2 

7 2 



+ 



ri + a 2 



r 2 ^ + a 2 cos 2 # 
With further simplifications, this becomes 

sin 2 fl 



\j3r + sin 9d9 — a sin 



;„2n ,71 2 



(24) 



ds 2 



t=0, r=r+ 



« + a 2 cos 2 9)d9 2 + 



[-Par+ sin 9d9 + (r 2 + a 2 )#] 2 . (25) 



ri + a 2 cos 2 



In the a — > limit, Eq. (25) gives precisely the expression we obtained for the boosted 
Schwarzschild metric. Let us now look at the 2-metric qab{A,B = 9,<p) for the apparent 
horizon component by component. 

\ 2 

7/3ar + sin 2 9 




ri + a 2 cos 2 9 
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ri + a 2 ) sin9 



r\ + a 2 cos 2 9 



r\ + a 2 ) sm9 
ri + a 2 cos 2 # 




For any 2x2 matrix of the form 



7/3ar + sin 2 9 



r\ + a 2 cos 2 6* 



(26) 



AB 



A 2 + B 2 BC 
BC C 2 



(27) 



the determinant is detHAB — AC. The 2-dimensional metric is of this form, so 



\J det (Jab) = ( r + + o 2 ) sin 9 



(28) 



Since 



Area = J a/ 'det (g AB )d9d(j) (29) 

we obtain an area of 47r fr 2 + a 2 ) as expected, identical to the unboosted horizon area. 

Going back to our boost in an arbitrary direction, we rewrite the 3-metric in Eq. (19) 
using the spheroidal coordinates of Eq. (21). After some algebra using a well known algebraic 
relation for the Kerr spacetime (2Mr + = rj. + a 2 ) to simplify, and setting f = r + in most 
places, we end up with a result surprisingly similar to Eq. (25) 

sin 2 fl 



ds\_ 



:0,r=r_)_ 



(r 2 + + a 2 cos 2 9) dW 



ri + a 2 cos 2 9 



[a^pdrw + (r\ + a 2 ) 



(30) 



Using the last one of Eqs. (21), we now expand the terms containing dr\\ and obtain the 
components of the 2-metric for the apparent horizon: 



(r 2 + a 2 cos 2 9) 

(3 2 a 2 sin 2 9 
-I — = 

ri + a 2 cos 2 9 



(r 2 + a 2 ) sin 2 9 
ri + a 2 cos 2 9 



a 2 sin 9 f3 cos 6* cos(0 — (j)^) — r + cos #3 sin 9 J , (31) 



+ a 2 — (3a sin 6^ sin 9 sin(</> — 0^) 



/^a/V 2 ,. + a 2 sin 2 9 
ri + a 2 cos 2 # 



/r + + a 2 sin 6*^ cos 9 cos (0 — 0^) — r + cos 6*^ sin 9 
^r\ + a 2 — (3a sin sin # sin(0 — 0^) 



(32) 
(33) 
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The /? — ^ limit of equations (31) through (32) yields the standard 2-metric of the Kerr 
spacetime given in Boyer-Lindquist coordinates. The a — > limit gives the standard 
Schwarzschild (spherical) 2-metric. The 8/3 = limit yields the metric of Eq. (26). Eqs. 
(31)-(33) show that g^B is again of the form of Eq. (27). Hence the square root of the 
determinant of the metric is 

\J det (<7ab) = {t\ + a 2 ) sin 9 — f3a\J r\ + a 2 sin 9p sin 2 # sin (0 — <f>p) . (34) 

The first term above is the familiar contribution from the unboosted Kerr metric. To 
determine the area, we integrate the square root of the determinant of the 2-metric over 
the angular variables of the spheroidal coordinate system. 

Area = j a/ det (gAB)d9d(f) 

(r 2 + + a 2 ) sin 9 — j3ayj r\ + a 2 sin 9p sin 2 9 sin (0 — (pp) 
= 4vr(r 2 +a 2 ) . (35) 

Above, the second term disappears because of the <f> integral. Our calculation shows that the 
area of the apparent horizon of a Kerr black hole remains invariant under arbitrary Lorentz 
boosts, as expected. 




IV. CONCLUSIONS 



Our goal was to show that the area of the apparent horizons of Kerr black holes remain 
invariant under a particular redefinition of the t = constant hypersurface (Lorentz boosts 
in arbitrary directions on the Kerr-Schild form). We introduced boost-parallel cylindrical 
coordinates. In this form, it is almost trivial to boost the spacetime metric. Once boosted, 
we looked at a t = constant hypersurface to determine the three dimensional spatial portion 
of the boosted metric, and projected down to the 2-metric of the apparent horizon. We gave 
examples and validation of using these results based on computational initial data, to obtain 
binding energy results for nonspinning black holes boosted together or apart. The binding 
energy curves with different (3 values closely overlap and agree in the Newtonian limit; 
the slight deviation of the binding energy from the 1/r form can plausibly be explained 
by nonlinear corrections to the physical separation corresponding to a given coordinate 
separation. This is a subject of a post-Newtonian study in progress. 
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For the general direction boosted Kerr case the 2-metric of the apparent horizon has 
non-zero off-diagonal terms; however, when integrated over the angular variables, these 
contribute zero to the area, leaving us with the same result as the undisturbed Kerr case, 
namely Area = An (r+ + a 2 ) . 

We performed all of our calculations in one particular type of slicing of the spacetime, 
i.e. t = constant. There are infinitely many slicings of static and stationary spacetimes, 
including ones with no apparent horizon at all [29]. But if the t = constant space contains 
an apparent horizon then the horizon area has the standard value. Black holes are very 
special objects indeed. 

Acknowledgments 

This work was supported by NSF grant PHY-0354842 and NASA grant NNG 04GL37G. 
Sarp Akcay would like to thank Cihan Akcay for his assistance with MATlab. 



[1] C. DeWitt, B.S. DeWitt, "Black Holes: Proceedings of the 23rd Les Houches Summer School", 

Gordon and Breach, New York (1973). 
[2] S.W. Hawking, G.F.R. Ellis, Large Scale Structure of Space-Time (Cambridge University 

Press, Cambridge, 1973). 
[3] Binary Black Hole Grand Alliance, Phys. Rev. Lett. 80, 2512-2516, (1998). 
[4] R.P. Kerr, A. Schild, "A New Class of Vacuum Solutions of the Einstein Field Equations", 

in Proceedings of the Galileo Galilei Centenary Meeting on General Relativity, Problems of 

Energy and Gravitational Waves, ed. G. Barbera (1965). 
[5] Mijan Firdous Huq, "Apparent Horizon Location in Numerical Spacetimes" Ph.D. Thesis, 

The University of Texas at Austin, (2004). 
[6] R. Matzner, M. F. Huq and D. Shoemaker, Phys. Rev. D59, 024015, (1999)[arXiv:gr- 

qc/9805023]. 

[7] M. F. Huq, M. Choptuik and R. A. Matzner, Phys. Rev. D66, 084024, (2002). [arXiv:gr- 
qc/0002076]. 

[8] E. Poisson, A Relativist's Toolkit (Cambridge University Press, Cambridge, 2004). 



17 



[9] C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravitation (W.H. Freeman, New York, 
1970). 

[10] M.P. Hobson, G. Efstathiou, and A.N. Lasenby, General Relativity (Cambridge University 

Press, Cambridge, 2006). 
[11] S. Chandrasekhar, The Mathematical Theory of Black Holes (Oxford University Press, Oxford, 

1983). 

[12] F. Pretorius, Phys. Rev. Lett. 95, 121101, (2005) [arXiv:gr-qc/0507014]. 

[13] J.G. Baker, J. Centrella, D. I. Choi, M. Koppitz and J. van Meter, Phys. Rev. Lett. 96, 111102, 

(2006) [arXiv:gr-qc/0511103]. 
[14] M. Campanelli, C. O. Lousto, P. Marronetti and Y. Zlochower, Phys. Rev. Lett. 96, 111101, 

(2006) [arXiv:gr-qc/0511048]. 

[15] J. A. Gonzalez, M. D. Hannam, U. Sperhake, B. Brngmann and S. Husa, arXiv:gr-qc/0702052, 

(2007) . 

[16] M. Campanelli, C. O. Lousto, Y. Zlochower and D. Merritt, Phys. Rev. Lett. 98, 231102, 

(2007) [arXiv:gr-qc/0702133]. 
[17] B. Brugmann, J. A. Gonzalez, M. Hannam, S. Husa and U. Sperhake, arXiv:0707.0135, (2007). 
[18] S. Brandt, B. Brugmann, Phys. Rev. Lett. 78, 3606, (1997). 
[19] R. A. Matzner, M. J. Huq, D. Shoemaker, Phys. Rev. D59, 024015, (1998). 
[20] R. A. Matzner, A. Nerozzi, P. Walter, "openGR", 2008, in preparation. 
[21] E. Bonning, P. Marronetti, D. Neilsen, R. A. Matzner, Phys. Rev. D68, 044019, (2003). 
[22] S. Hawley, M. Vitalo, R. A. Matzner, arXiv:gr-qc/0604100, (2006). 

[23] R. Genzel, R. Schodel, T. Ott, A. Eckart, T. Alexander, F. Lacombe, D. Rouan and B. 
Aschenbach, Nature 425, 934, (2003). 

[24] R. Shafee, J. E. McClintock, R. Narayan, S. W. Davis, L. Li, R. A. Remillard, The Astro- 
physical Journal 636, L113, (2006). 

[25] A. Broderick, V. L. Fish, S. S. Doeleman, A. Loeb, arXiv.org: 0809.4490, (2008). 

[26] R. Matzner, "Area of the Apparent Horizon for a Schwarzschild Black Hole is Invariant under 
Lorentz Boosts in the z-direction." unpublished, (1996). 

[27] R.E. Wald, General Relativity (The University of Chicago Press, Chicago, 1984). 

[28] R. Arnowitt, S. Deser, and C. Misner, Gravitation, an Introduction to Current Research (Wi- 
ley, New York 1962). 

18 



[29] R.E. Wald, V. Iyer, Phys. Rev. D44, 3719, (1991). 



19 



